Phase transitions, entanglement and quantum noise interferometry in cold atoms 
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We show that entanglement monotones can characterize the pronounced enhancement of entanglement at 
a quantum phase transition if they are sensitive to long-range high order correlations. These monotones are 
found to develop a sharp peak at the critical point and to exhibit universal scaling. We demonstrate that similar 
features are shared by noise correlations and verify that these experimentally accessible quantities indeed encode 
entanglement information and probe separability. 



Quantum entanglement represents one of the most fascinat- 
ing features of quantum theory and has emerged as an im- 
portant resource in quantum information science [1J. As en- 
tanglement represents a unique form of correlations that do 
not occur in classical systems, the investigation of connec- 
tions between entanglement and quantum phase transitions is 
an emerging field of research which promises to complement 
the understanding of critical phenomena in condensed matter 
physics and quantum field theory. 

Iconic examples of exactly solvable models exhibiting 
quantum phase transitions, such as the Ising-like spin chain, 
are becoming paradigms for investigating the suitability of 
proposed entanglement measures, e.g. the concurrence and 
entanglement entropy [2 J, to quantify the entangling resources 
of a system close to its quantum critical point. Among the var- 
ious entanglement measures the concurrence has [2| gained 
particular attention in view of its universal scaling OHH close 
to the phase transition. However, the regular bipartite concur- 
rence has the drawback that even at the critical point where 
spin-spin correlations extend over a long range (as the cor- 
relation length is diverging for an infinite system) only the 
next and next-to nearest neighbor concurrences are non zero. 
Moreover, it is not the concurrence that peaks at the transition 
but its first derivative. 

The expected enhancement of entanglement at the phase 
transition has been observed with the help of generalizations 
of the bipartite concurrence or the entanglement entropy 
St between a block of £ consecutive spins and the rest of the 
chain [6, 7]. However, there is no systematic understanding of 
why some entanglement measures capture this enhancement 
whereas other fail to do so. In the present letter we provide a 
general explanation for the different abilities of various entan- 
glement measures to describe quantum phase transitions. 

With the help of a sequence of entanglement measures we 
illustrate the importance of higher order correlations in en- 
coding long range correlations near quantum phase transition. 
For this purpose we characterize the entanglement properties 
of a chain of interacting spins with different entanglement 
monotones and study how correlations are established over 
different ranges. We see that there is a qualitatively different 



behavior of long-range correlations as opposed to their short- 
range analogues, as well as between two- and many-body en- 
tanglement. 

Whereas entanglement measures have become common 
tools for theoretical investigations, they are significantly less 
popular in the experimental community since they are typi- 
cally extremely arduous if not practically impossible to mea- 
sure. For the verification of theoretical findings and for the in- 
vestigation of theoretically intractable systems it is therefore 
necessary to have efficient experimental means to probe en- 
tanglement properties. Among the various physical systems 
exhibiting quantum phase transitions, cold atoms loaded in 
optical lattices offer considerable advantages over more tra- 
ditional set-ups since they allow to engineer various Hamil- 
tonians with high precision and to include or remove noise 
and randomness in a controlled way. In these systems many- 
body correlations are extractable from absorption images of 
the atomic cloud after its release from the trap J8j |9] [lOl . 
However such images typically do not provide sufficient in- 
formation to recover a proper entanglement measure. Here we 
give a general prescription of how to extract the pronounced 
enhancement of entanglement at a phase transition from the 
limited accessible experimental data. 

We consider the ground state of the one dimensional spin- 
1/2 anisotropic XY model in a transversal magnetic field 
characterized by the Hamiltonian 
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where J is the nearest neighbor coupling constant, h is a trans- 
verse magnetic field, 7 is the anisotropy parameter, < 7 < 
1, er" are Pauli matrices (a = x,y,z) and periodic bound- 
ary conditions are assumed throughout. For < 7 < 1, H 
belongs to the Ising universality class and exhibits a quan- 
tum phase transition from a paramagnetic to a ferromagnetic 
phase when A = h/J takes its critical value A c = 1 ifTTI . The 
XY model is exactly solvable, and any correlation function 
can be expressed in terms of Toplitz-like determinants after 
a Jordan- Wigner transformation that maps spin operator into 
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fermionic operators |[T2l \l3l . The evaluation of all four-point 
correlations is however more cumbersome than the commonly 
used prescription for computing two-point functions, and de- 
tails about the procedure will be provided elsewhere[21 1. 

A severe problem for the investigations of entanglement in 
many-body systems is that such entanglement often cannot 
be discovered in terms of few -body correlations. Thus, a few- 
body state that is obtained via a partial trace over a many -body 
state often shows only weak or no entanglement even if the 
underlying many -body state is highly entangled EJHJ. There- 
fore, we will not consider entanglement properties of reduced 
few-body states, but only investigate properties of the ground 
state of the entire many-body system. Instead of focusing on 
quantifying the entanglement remaining in a subsystem of a 
few spins after tracing over the others, we focus on the entan- 
glement shared between these spins and the rest. The latter 
does encode information about the entanglement of the pure 
many-body state. 

There is an enormous zoo of entanglement monotones and 
measures to quantify the entanglement of many-body sys- 
tems lfl4l . and there is no generally accepted choice that el- 
evates one above the rest. For bipartite systems, both en- 
tanglement entropy and concurrence play an important role, 
since they are the first monotones that could be evaluated by 
purely algebraic means for mixed states lfT31 [161 . Entangle- 
ment entropy is defined in terms of the degree of mixing of 
the reduced density matrix g r of either of the two subsystems: 
S(^) — — Trg r log g r , and the most frequently used defini- 
tion of concurrence is c(*f?) = 1(^*1^ <g> <7 y |'I')|, where a y is 
the second Pauli matrix, ("f* | is the complex conjugate of 
that is the transpose of |\&), with the conjugation/transposition 
performed in the eigenbasis of the third Pauli matrix a z . For 
a system of two spin- 1/2 objects concurrence and entangle- 
ment entropy are equivalent as one is a bijective function of 
the other. However, for the case of multipartite systems that 
we consider here, this is no longer the case; in particular, there 
is no unique generalization to many -body systems for either of 
these measures, and our choices are taken in accordance with 
the currently investigated situation without any claim to be a 
canonical choice. 

In the case of entanglement entropy we will divide the spin 
chain into two blocks of not necessarily consecutive spins and 
consider the entanglement entropy between those blocks. In 
the case of concurrence we will use generalizations based on 
the popular redefinition c(\E') = -\/2(l — Tig^) [ 1 7 1 of con- 
currence: we will focus on the discrete set of generalized tan- 
gles (squared concurrence) lfl8l : 
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where the q v are reduced density matrices of 1 to k spins, and 
the summation is performed over all reduced density matrices 
of up to k spins. Formally, this means that contains all 1 to 
k touples of pairwise different elements labeling the individ- 
ual sites. The constants Dj~ — J2i=o TVj=o(N~j) are chosen 
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FIG. 1: (color online) Various multipartite variations of the tangle 
(squared concurrence) Ti long-dashed (pink), T2 dashed (yellow), 
T3 solid(red), and T4 dotted (blue) as function of A for 7 = 1 (right 
panel) and 7 = 0.5 (left panel). Ti, and T2 do not show a distinct 
peak around the phase transition. T k are obtained by numerical com- 
putation of Toplitz determinants for a system of size N = 175. 



such that fully separable states yield a vanishing value of T k . 
The advantage of this choice of entanglement monotone as 
compared to the various alternative options is that the present 
objects can easily be evaluated in terms of the four point cor- 
relation functions that can be expressed in terms of Toplitz 
determinants. T\ describes the average entanglement between 
an individual spin and the rest of the chain. T2 includes in 
addition the entanglement contribution between pairs of spins 
and the residual system, and T 3 , and T 4 contain ever higher 
order correlations. This hierarchy extends to the multiparticle 
concurrence |fl9l that provides a quantification for the entire 
entanglement content of a multipartite system. However, for 
practical reasons we will restrict the present investigation to 
T1-T4 since the evaluation of 5-point quantities becomes im- 
practical for large systems. 

A clear difference between the behavior of bipartite and 
multipartite entanglement can be seen already going from T± 
to T4 as shown in Fig.[T[ where the different multi-partite tan- 
gles are plotted as function of A for 7 = 1, and 7 = 0.5. 
In both cases neither T\ nor T2 show a distinct peak around 
the phase transition, but their behavior rather resembles the 
step-like behavior of the bipartite concurrence [3]. But start- 
ing from T 3 for the case of 7 = 1, and T 4 for 7 = 0.5, 
there is a clear peak arising around the critical point. As it 
can be seen for 7 = 1 the peak becomes more pronounced 
with increased order of correlation in T k . This shows that the 
long-range correlations that are established around a quantum 
phase transition are displayed in terms of multipartite entan- 
glement, whereas the bipartite entanglement between blocks 
of consecutive spins do not display this behavior. 

In order to quantify the growth of multipartite entanglement 
in T4 at the critical point we study the scaling with system size 
of its peak height and position. The right inset of Fig.|2]shows 
the growth of the peak with increasing number of spins. Since 
T4 is a bounded quantity, it can not diverge even in the ther- 
modynamic limit. Nevertheless, the increase follows a log- 
arithmic fit (solid line). The left inset of Fig. [2] shows that, 
with increasing system size, the position of the peak, A m , 
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FIG. 2: The main plot shows the universality of V = (T 4 — T4 ) _1 
around the phase transition: T(A) — T(A m ) is plotted for various 
system sizes, from TV = 75 to 175, as function of TV(A — A m ), and 
the same behavior is found for all system sizes that are displayed with 
different symbols. The insets show the growth (right) and position 
(left) of the peak with increasing system size. 



gets shifted towards the expected value in the thermodynamic 
limit, X m = 1, however, according to the logarithmic fit that 
follows the data, it does not reach A m = 1. Instead it over- 
shoots a bit towards the value of 0.97. We do not attribute 
too much meaning to this observation, but rather expect that 
more data for larger systems would result in a slightly more 
accurate extrapolation. 

Considering the well known universal behavior of Ising- 
like models at a phase transition, it is central to check whether 
our presently utilized quantification of entanglement reveals 
universality. Since by definition T4 is always finite a more 
appropriate quantity to study universal scaling is T = (T 4 — 
T|) _1 since it exhibits a logarithmic divergence at the crit- 
ical point. Here T 4 is the maximum value reached by T 4 
in the thermodynamic limit at A m = 1, which we calculate 
to be 7r/2. Again T| does not agree exactly with the value 
we get from our finite size scaling which yields 1.77 how- 
ever as we mention before we attribute it to the size limi- 
tation we have in calculating four point correlations and ex- 
pect a more accurate fit for larger systems. Fig. [2] shows that 
r(A) - r(A m ) = - A m )] exhibits universal be- 

havior with the scaling exponent v = 1 characteristic for the 
Ising model [11 J: When we plot this quantity as function of 
TV(A — A m ) for different system sizes, there is an almost per- 
fect coincidence of the various data points into a single curve. 
Furthermore, the shape of the universal curve closely resem- 
bles the universal curve underlying concurrence Q. 

In order to show that it is exactly the long range correla- 
tions that grow, as opposed the entanglement between close 
neighbors, and that the above findings are not a particularity 
of the employed tangles Ti, we present a similar analysis as 
above also for a different, well established entanglement mea- 
sure, the entanglement entropy. In contrast to previous stud- 
ies (6l |7] which considered bipartition of the spin chain into £ 
consecutive spins and the rest of the chain here we study the 



entanglement entropy S^L) extracted by bipartition of the 
spin-chain into four spins that are separated by equal distance 
L, (say spin 1, 1 + L, 1 + 2L, and 1 + 3L) and the rest of the 
system. This quantity is shown in Fig. [3] Note that similarly 
to T 4 , also 5 4 (L) remains finite in the thermodynamic limit, 
Si(L) < 4. While there is only a tiny enhancement of the 
entanglement entropy around the phase transition for L = 1, 
a peak grows with increasing separation L of the spins, and 
takes its maximum for L ~ TV/4, i.e. in the case of maxi- 
mum separation. This behavior again gives evidence of the 
long-range character of the correlations. 
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FIG. 3: Entanglement entropy, 5*4(1/), as function of A for differ- 
ent spacings between the four spins. 5*4(1/) grows monotonously 
with increasing separation and reaches its maximum value when 
L = TV/4. Here TV = 201. 

Having defined proper entanglement quantifiers that cap- 
ture the properties of a quantum phase transition we will now 
address the issue of how the latter can be linked to experimen- 
tally accessible observables. Here we will focus our discus- 
sion on cold atomic systems in view of their appeal as ideal 
quantum simulators of iconic condensed matter Hamiltoni- 
ans. In these systems absorption images taken after releas- 
ing the atoms from the trap are the commonly used diagnos- 
tic tools: the average density profile after time of flight maps 
to the quasi-momentum distribution n(q) of the atoms at the 
release time, and the density-density correlations known as 
noise correlations [8] yield the momentum-momentum corre- 
lations A(q 1 ,q 2 ) . These functions are just Fourier transforms 
of two and four-point correlations of the atomic field creation 
and annihilation operators at the various lattice sites, a n , : 



n(q) 
A(gi,<? 2 ) 



1 



N ^ 

n,m 



q(n-m)/tf^\ (3) 



\CL n Cl 7n 

{(n qi " (n qi ))(n qi - (n qi ))}) . (4) 



The noise correlations A(<ft, q 2 ) contain the same non-local 
correlations as those included in T 4 and S* 4 (i). In contrast 
to T 4 and S4, A (qi, q 2 ) is not a proper entanglement mono- 
tone, but it is accessible to measurement. Therefore, the nat- 
ural question that we address below is: how much informa- 
tion about the enhancement of many-body entanglement at the 
critical point can be inferred from noise correlations. 
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In view of the well known mapping between hard-core 
bosons (HCB) and spin- 1/2 operators iTOl . noise correlations 
for the spin chains can be defined by using the following re- 
lations between the spins and atomic operators: er+ = at, 

aj = hj and cr| = a/jCtj. In the following we use this map- 
ping to study the behavior of A(0, 0) and n(0) as A is varied 
across the critical point. Our aim is to use this exactly solvable 
system as a benchmark of the behavior of these correlations in 
more general systems, which are harder to deal theoretically 
but on the other hand where noise correlations can be accessed 
experimentally. Whereas their quantitative properties will de- 
pend on the details of the underlying Hamiltonian, we believe 
their general features will be generic due to the inherent uni- 
versal scaling of a quantum phase transition. 

Fig. |4]shows that A(0, 0) captures the enhancement of en- 
tanglement at the quantum critical point, whereas n(0) fails 
to do so. While n(0) exhibits a step-like functional depen- 
dence of A (similar to T^), the noise autocorrelation function 
A(0, 0) is sharply peaked around the phase transition (similar 
to T4). A (0,0) contains the higher order non-local correla- 
tions which are the key to properly extract the enhancement 
of entanglement at the quantum phase transition. 
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FIG. 4: (color online) n(0) and A(0, 0) as a function of A for N = 
95 spins and two different anisotropies: The dot-dashed blue (red 
dotted) curves and solid black (dashed green) correspond to n(0) and 
A(0, 0) for 7 = 1(7 = 0.5). The horizontal line at 12 is the maximal 
value that A(0, 0) can take for a separable state and the peaks of 
A(0, 0) exceed this threshold. The inset shows the numerical scaling 
of the A(0,0) peak vs N which is in agreement with the analytic 
prediction of a powerlaw with exponent 3/2. 

Our calculations for noise correlations also verify that they 
can be related to entanglement because they probe separabil- 
ity. For any separable state |$ s ) of a spin-1/2 system (i.e. a 
state that can be expressed in terms of one-body states |<^) as 
\$s) = |<Ai) ® \h) ® • • • ® |<MX yields a value of A(0, 0) 
that is smaller than ED 

Amax($ s ) = g (1 + N) (5) 

That is, if the obtained value of A(0,0) exceeds this thresh- 
old then the underlying state it entangled. In Fig. [4] we indi- 
cate Amax(Ps) for N = 95 . The peak of A(0, 0) clearly 



exceeds the separability threshold and the excess becomes 
even more pronounced with increasing system size. The rea- 
son is the following: While A(0,0) grows algebraically with 
the system size as iV 3 / 2 , the separability threshold of Eq. |5]) 
grows linearly with N, so that the excess of the peak over this 
threshold increases as y/N. The exponent 3/2 in the noise 
correlation peak height follows from the universal scaling of 
(ct x ) ~ (A — l) 1 / 8 which yields a divergence of the noise cor- 
relation peak as (A — l) -3 / 2 in the thermodynamic limit and 
to a N 3 / 2 scaling for finite samples. 

Even though the above investigations were focused on the 
anisotropic XY model Hamiltonian, we emphasize that Eq. 
(5) provides a general criterion for arbitrary spin systems and 
consequently corresponds to a useful benchmark to study en- 
tanglement properties in more general systems such disor- 
dered spin or spins at finite temperature. Additionally, pre- 
liminary calculations done in in comparatively small ID soft 
core bosons undergoing a superfluid to Mott insulator transi- 
tion lfT2H20]| also indicate a similar peaked behavior of A close 
to the critical point. This suggests that the growth of A can 
be a generic signature of a quantum phase transition and that 
noise correlations are suitable observables to experimentally 
verify the enhancement of entanglement. 
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